Lower Memory Oblivious (Tensor) Subspace Embeddings with Fewer Random Bits: Modewise Methods for Least Squares
نویسندگان
چکیده
In this paper new general modewise Johnson--Lindenstrauss (JL) subspace embeddings are proposed that can be both generated much faster and stored more easily than traditional JL when working with extremely large vectors and/or tensors. Corresponding embedding results then proven for two different types of low-dimensional (tensor) subspaces. The first these produces improved space complexity bounds rank-$r$ tensors whose CP decompositions contained in the span a fixed (but unknown) set $r$ rank-$1$ basis vector setting result yields very near-optimal oblivious constructions require fewer random bits to generate standard subspaces $\mathbb{C}^N$ spanned by special Kronecker structure. second herein provides fast arbitrary $r$-dimensional $\mathcal{S} \subset \mathbb{C}^N$ which also (and so easier store, i.e., less space) methods order achieve small $\epsilon$-distortions. These work (i) effectively folding any given $\mathcal{S}$ into (not necessarily low-rank) tensor, (ii) resulting tensor $\mathbb{C}^m$ $m \leq C r \log^c(N) / \epsilon^2$. Applications related compression compressed least squares solution considered, including those used fitting low-rank decompositions, shown well numerically settings.
منابع مشابه
Lower Bounds for Oblivious Subspace Embeddings
An oblivious subspace embedding (OSE) for some ε, δ ∈ (0, 1/3) and d ≤ m ≤ n is a distribution D over Rm×n such that for any linear subspace W ⊂ Rn of dimension d, P Π∼D (∀x ∈W, (1− ε)‖x‖2 ≤ ‖Πx‖2 ≤ (1 + ε)‖x‖2) ≥ 1− δ. We prove that any OSE with δ < 1/3 must have m = Ω((d + log(1/δ))/ε2), which is optimal. Furthermore, if every Π in the support of D is sparse, having at most s non-zero entries...
متن کاملPrimality Testing with Fewer Random Bits 1
In the usual formulations of the Miller-Rabin and Solovay-Strassen primality testing algorithms, to test a number n for primality, the algorithm chooses \candidates" x 1 ; x 2 ; : : :; x k uniformly and independently at random from Z n , and tests if any are a \witness" to the compositeness of n. For either algorithm, the probability that it errs is at most 2 ?k. In this paper, we study the err...
متن کاملTight Bounds for $\ell_p$ Oblivious Subspace Embeddings
An lp oblivious subspace embedding is a distribution over r × n matrices Π such that for any fixed n× d matrix A, Pr Π [for all x, ‖Ax‖p ≤ ‖ΠAx‖p ≤ κ‖Ax‖p] ≥ 9/10, where r is the dimension of the embedding, κ is the distortion of the embedding, and for an n-dimensional vector y, ‖y‖p = ( ∑n i=1 |yi|) 1/p is the lp-norm. Another important property is the sparsity of Π, that is, the maximum numbe...
متن کاملInner-Iteration Krylov Subspace Methods for Least Squares Problems
Stationary inner iterations in combination with Krylov subspace methods are proposed for least squares problems. The inner iterations are efficient in terms of computational work and memory, and serve as powerful preconditioners also for ill-conditioned and rank-deficient least squares problems. Theoretical justifications for using the inner iterations as preconditioners are presented. Numerica...
متن کاملParallel Tensor Methods for Nonlinear Equations and Nonlinear Least Squares
We describe the design and computational performance of parallel row-oriented tensor algorithms for the solution of dense systems of nonlinear equations and nonlinear least squares problems on a distributed-memory MIMD multiprocessor. Tensor methods are general purpose methods that base each iteration upon a quadratic model of the nonlinear function, rather than the standard linear model, where...
متن کاملذخیره در منابع من
با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید
ژورنال
عنوان ژورنال: SIAM Journal on Matrix Analysis and Applications
سال: 2021
ISSN: ['1095-7162', '0895-4798']
DOI: https://doi.org/10.1137/19m1308116